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Abstract 

The mechanisms of energy and spin transfer in quantum dot pairs coupled via the Coulomb 
interaction are studied. Exciton transfer can be resonant or phonon-assisted. In both cases, the 
transfer rates strongly depend on the resonance conditions. The spin selection rules in the transfer 
process come from the exchange and spin-orbit interactions. The character of energy dissipation 
in spin transfer is different than that in the traditional spin currents. The spin-dependent photon 
cross-correlation functions reflect the exciton transfer process. In addition, a mathematical method 
to calculate Forster transfer in crystalline nanostructures beyond the dipole-dipole approximation 
is described. 

PACS numbers: 78.67.Hc, 72.25.Fe, 73.21.La 
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A new field of research, spintronics, studies the principles of manipulation of the spin 
degree of freedom in solids and molecules^, whereas traditional electronics utilizes electric 
charges. Spintronics is also closely related to quantum information science since the spin 
is an important element of quantum computing. In most cases, transport of spins in solids 
and molecular systems comes from transfer or tunneling of charged electrons and, therefore, 
is accompanied by electric currents. In electronic materials, the electric interactions are 
often much stronger than the spin-related ones and, therefore, when manipulating charged 
particles with spins, usual electronics has often obvious advantages compared to spintron- 
ics. However, spin may have advantages over charge. In contrast to charge or mass, the 
angular momentum or spin can be transferred without tunneling or ballistic transport. One 
particular mechanism of spin transfer without tunneling can occur in optically excited semi- 
conductor quantum dots (QDs); spin-polarized excitons can be transferred between QDs 
via the long-range noncontact Coulomb interaction 2 ". It is important to note that Coulomb 
(Forster) transfer of spin in QDs becomes possible due to the strong spin-dependent inter- 
actions in semiconductors, such as spin-orbit and exchange couplings 2 . 

Here we study theoretically the microscopic mechanisms of spin-dependent Forster trans- 
fer in a molecule composed of two self-assembled QDs (figs. 1 1121) . In the typical scheme of 
Forster transfer—, an optically excited exciton in QD1 ("donor") becomes transferred to QD2 
("acceptor") via the Coulomb interaction (fig. |2k). The traditional methods to observe this 
transfer are time-resolved photoluminescence (PL) spectroscopy^ and photon correlations 5 . 
In the case of resonant transfer in self-assembled QDs, the spin selection rules are determined 
by the electron-hole exchange interaction in an exciton and by the spin-orbit interaction in 
the valence band. In a symmetric QD molecule, transfer occurs with conservation of the ex- 
citon spin configuration, whereas, in QD molecules with broken symmetry, the exciton spin 
becomes partially lost in the transfer process. The transfer rates exhibit a strong dependence 
from the exciton energy difference in a QD pair, AE = E exc ^oti — E exc ^ot2- In the resonant 
regime AE ~ 0, exciton and spin transfer is fast. In the nonresonant regime, transfer can 
be assisted by acoustic phonons. Again, it strongly depends on AE. In contrast to the 
previous paper on spin transfer—, we include here the electron-hole exchange interaction and 
consider the fine structure of exciton states. In addition, we show that the dipole-dipole 
approximation is not reliable for the typical inter-dot distances in experimental structures 
and describe a method to compute Coulomb matrix elements beyond the dipole approach. 
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Our method is valid when R ^> aiattice, where R is the inter-dot distance and aiattice is the 
lattice period. In addition, we note that the Forster transfer mechanism considered in this 
paper has the electrostatic, near-field nature; this is in contrast to the recent paper on the 
radiative coupling between QDs 6 . 

Experimentally, Forster transfer of excitons can be studied using time-resolve 
photoluminescence^ or the photon cross-correlation metho d 5 ^ 8 . Experiments on energy 
transfer in nano-structures were done with colloidal nano-crystals^ and recently with self- 
assembled InAs QDs^. It is important to note that QDs can also be coupled via tunneling^. 
However, the tunneling amplitude decreases exponentially with increasing the distance be- 
tween QDs and with the hight of the barrier between QDs. At the same time, the Forster 
transfer rate demonstrates a power law (oc R~ 6 ) and is independent of the hight of the poten- 
tial barrier between QDs. To avoid tunneling in a QD molecule, one can grow AlGaAs barrier 
between the QDs or one can use QDs with stronger confinement. Another important factor 
is that the resonance conditions for Forster and tunneling transfers are different. Therefore, 
it seems to be possible, by using self-organization growth, to design self-assembled QDs 
with predominantly Forster transfer. The recent experiment^ indicates that indeed such 
pairs of QDs with Coulomb-induced coupling can be designed using the self-organization 
growth method. Another recent experimental paper (ref.— ) describes the spin-response of 
colloidal QDs bridged with bio-molecules which may assist direct tunnel transport between 
nano-crystal o 1Ql11 . 

The spin current of mobile polarized electrons is always accompanied by the electric 
charge flow. This brings back the old issue of energy dissipation in electronic devices, 
given by the Joule heat, Q = jE (here j is the electric current and E the driving electric 
field). In the case of Coulomb-induced transfer between QDs (figs. |2t and|§J), the energy 
dissipation has a different character and comes from phonon-assisted relaxation. The energy 
loss in this process is equal to the energy level difference in the donor and acceptor QDs: 
AE = E exCt doti — E exc ^ot2 = TuJph- In such processes, the energy AE turns into the phonon 
energy hu p h. If a pair of QDs is resonant (AE ~ 0), the coupling between QDs can become 
coheren t 12 i 13 i 14 ; such a process of spin transfer resembles Rabi oscillations between QDs 
(fig. Et). Coherent spin transfer can occur without dissipation. 

The paper is organized as follows: the section 1 describes the model of the QD system, 
the section 2 includes the results on spin-dependent transfer rates, the section 3 is devoted 
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FIG. 1: (a) Schematics of the systems of two QDs. (b) Geometry of a pair of self-assembled QDs 
and the corresponding crystallographic axes, (c) Geometry of two QDs with broken symmetry; the 
vector d describes the shift. 

to the phonon-assisted transfer, and the sections 4 and 5 are about the photon correlation 
functions and transfer under strongly resonant conditions. 



I. MODEL 

We now consider a model of a pair of self-assembled QDs^^ ? . Our model incorporates 
two oblate asymmetric QDs (figs, and b); the vertical, z-size of QDs is assumed to be 
much less than the lateral ones. To model the lateral motion of electrons and holes, we use 
the harmonic functions^ with the characteristic lengthes, l e (h),x — y fr/wx^m e (h) and $® = 
\J h/ujy^m^h), where are the harmonic-oscillator frequencies for electrons (e) and holes 
(h), m e (h) are the effective masses of particles, and (x,y) are the in-plane coordinates. The 
complete envelope functions used below have a form: = A e ^e~ x2 ^ 2l ^ x ~ y2 ^ 2l ^ysin[ii(z — 
Zk)/L] and <ph,k = ^h,k€~ x ^ 2lh < x ~ y ^ 2lh >ysin\jr{z — Zk)/L], where i — x, y, k — 1, 2 is the dot 
number, Zk is the ^-coordinate of the center of fc-dot, Afc )e (M are the normalizing coefficients, 
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FIG. 2: Schematics of transfer processes between two QDs. (a) Phonon-assisted mechanism, (b) 
Resonant coherent transfer. 

and L is the "vertical" size of QDs. In the following, we will use the system of coordinates 
(x, y) (fig. [TJd) which corresponds to the typical orientation of elongated self-assembled QDs 
grown on the [001] surface^. The lowest excitonic states in QDs responsible for PL are 
composed of heavy holes and electrons and correspond to the s-like envelope wave function. 
Taking into account the only heavy-hole wave functions, we can write the spin Hamiltonian 
of exciton of an individual QD in the following forr^ ' | s 



#e-I = a zjzS z + b rf?Sii (!) 

i=x,y,z 

where s» is the electron spin matrixes and Sj are the 2x2 angular-momentum operators of 
heavy holes. The exchange parameters a z and hi depend on a particular QD. By using the 
operator (fljl. we find the exciton wave functions and their energies. The bright excitons are 
composed of the states with J tot = j z + s z = ±1: 



I + -• — - > 4-1 — -• 4-- > 

b _ \ ~ 2> 2 ~\ 2' ~ 2 



I + I._3 > _| _ l. + 3 > 

^= ' + 2 ' 2 J 2 ' +2 (2) 



where we used the notation \s z ;j z >. The corresponding energies e b x = — (|a 2 4- j^b z ) + 
§(frz — by) and e b y = —{\a z 4- j^b z ) — ~(b x — b y ). The lower indexes x, y reflect the character of 
spin orientation in an exciton and the optical selection rules. In the PL process, the excitons 
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ijj x and ij)y create photons with linear x and y polarizations, respectively. The dark excitons 
are composed of the states with J tot = j z + s z = ±2: 





(3) 



V2 



Their energies e b x = -(\a z + ^b z ) + §(6^ + b y ) and e b y = -{\a z + f|fe 2 ) - \{b x + b y ). Typically, 



the two lowest states in the exciton spectrum are dark whereas the two upper ones are bright 
(fig. 121). In this model, the normal magnetic field does not lead to mixing between dark and 
bright states, inducing an additional splitting in the pairs of states. In the limit B —>■ oo, 
the wave functions become the states in which the angular momentum is a good quantum 
number. 

The exchange spin-dependent interaction in excitons and the dark-bright energy splitting 
are quite strong, about 0.5 meV^. The other types of interaction in crystalline QDs, inter- 
dot Coulomb and electron-phonon couplings, can be weaker and we are going to involve 
them as perturbation. We note that the intra-dot Coulomb interaction is quite strong, but 
it does not lead to the inter-dot exciton transfer; it mostly shifts down the exciton energies. 
Then, the perturbation Hamiltonian is 



where Ucoui is the inter-dot Coulomb interaction. The operators H e _ p h and Hh- P h represent 
the interaction between acoustic phonons and particles. 

II. SPIN-DEPENDENT COULOMB MATRIX ELEMENTS 

First we compute the inter-dot Coulomb matrix elements. The complete set of electron- 
hole wave functions includes 8 states: 



Hp er turb — Ucoul{r\, T 2 ) + H e ^ p h + Hf^ph, 



(4) 



± 



(5) 
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where k = 1, 2 in the QD index. For the one-exciton states, we have \s' z ,j' z ; 1 > |0, 2 > and 
|0; 1 > 2 >; here |0; k > denotes the state of the fc-dot without an exciton. Then we 

write the inter-dot Coulomb matrix elements as: 

< 0;1| < s 2 ,j2;2\U C oui\si,ji; 1 > |0,2 > . (6) 

In most papers, the matrix elements (jBJ) are calculated within the dipole-dipole approxi- 
mation which is valid in the limit R -C Idoti where R is the inter-dot distance and I dot is a 
characteristic size of QDs {Idot ~ l e (h),x(y))- Now we are going to use a method beyond the 
dipole-dipole approximation. Namely, we are going to use the quantity aiattice/R as a small 
parameter, where aiattice is the crystal- lattice period. Since, aiauice *C Idot, our approxima- 
tion is much better compared with the standard dipole-dipole approximation, R ^> Idot- To 
evaluate the matrix element (JHJ), it is convenient to return to the pure electron representation 
and to consider two electrons in each QD explicitly: 

M\ Su j 1 -dotl>^\s 2 ,j 2 ;dot2> =< 0; 1| < S2,j2',2\Ucoul\Sl,jl', 1 > |0, 2 > = 
J ^^i^ 2 ^2[*i l! l(ei)*-i 1) l(ei)* S2) 2(6)% 2) 2(e2)]* 

(7) 

here £1 = (1*1,01) an d = (^,0^) are the spatial and spin coordinates of two electrons 
in QD1; £2 = (r 2 , 02) and £ 2 = (r^J ^) are similar coordinates for QD2. su2) are the z- 
components of spin in the conduction band of QD1(2); ji(2) are the z-components of angular 
momentum of the valence band-electrons. ^j kl k is the electron state in the valence band 
of the k-QD {k = 1,2, = ±3/2). The function ^ Skt k is the state for the conduction 
band (s^ = ±1/2). Within the envelope-function approach ^j^iO = Xj( r , °")0/i,j( r ) an d 
^s,k(0 — Xs( r , cr )4 ) e,i( r ), where Xj ari d Xs ar e the Bloch functions in the valence band 
(heavy- hole) and in the conduction band, respectively; <ph(e) ( r ) is the hole (electron) envelope 
function. The Forster process involves one electron in QD1 with coordinate r% and one 
electron in QD2 with coordinate r 2 (fig. EJ) and therefore we can integrate over rj^. For 
smooth envelope functions and long-range Coulomb potential, we can rewrite the integral 
(J2J) as follows: 
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FIG. 3: Electron configurations for the initial and final states of the transfer processes with (a) 
and without (b) conservation of the exciton angular momentum. The process (b) becomes possible 
in QD pairs with broken symmetry (d ^ 0). 



<Pe,l{R ai )<Ph,l{R ai )<Pe,2 (Ra 2 )<f>h,2 {Ra 2 ) 

* I I dAeidAe 2 [x jl (ei)*X S2 (6)*^c ou K^+Ai2 ai ,i2 a2 + Ai2 a J Xsi (6)x J - 2 (e2)], 

(8) 

where the summation is performed over all the unit cells in both QDs; a.\i2) are the unit 
cell indexes; Q ak and R ak are the unit cell volumes and unit cell coordinates, respectively 
(k = 1,2). A£ fc = (AR ah , o"fc), where AR ak is the spatial coordinate relative to the center of 
the a-cell and is the spin coordinate. Assuming ai a ttice/ldat ~ ^lattice I R < 1, we expand 
the Coulomb potential in terms of AR ak and take into account the leading term responsible 
for the Forster transfer: 
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tt ( r> i AD D i AD ^ e A-R^A-R^j 3 [( Ai? ai jQ , 2 ) ( A.R Q2 i? Ql jQ , 2 )] / 1 R ai )Q2 1 

tVCou/l-Kai + ^CaD /t a2 + Ait, 



" 2; e IE I 3 

c I ± L a\ ,02 I 



(9) 



where R ai , a2 = Ron ~ Ra 2 - The Coulomb potential was taken in the usual form: Uc ou i = 
e 2 /(e\ri — r 2 |). Changing the summation in eq. (JHJ) to integration, we obtain 

e 2 /" 

^|si,ii;dotl>-+|s 2 ,i 2 ;dot2> = ~ dRidR 2 (p e ^(Ri)(() h ^(Ri)(() e 2(R 2 )(j) h 2(R 2 ) 
[ [ r!A£ rl\£ \ (KfiY ^^ ^Rl^R2-^RlRl2)(AR2Rl2)]/\Rl2\ 2 . Afi . . Afi „ 

Jn, Jn 2 \ n i2\ 



(10) 



where R\ 2 — R\ — R2 and A£& = (A/4, <Jk). In the next step, we derive the matrix elements 
for the heavy-hole excitons using the Bloch functions, X1/2 = \S >|, X-1/2 — >|, 
X+3/2 = (\X > +i\Y >) T /y/2, X-3/2 = (\X > -i\Y >) | /y/2, 



Mq — M| Sl=+ i/2 ) j 1 =-3/2;dotl>^|s 2 =+l/2,j2=-3/2;(iot2> 
= M| Sl= _x/2J 1 =+3/2;dotl>-+|s 2 =-l/2Ja=+3/2;rfot2> 

1 _ 3 (Xi-.y 2 ) 2 +(yi-y 2 ) 2 
= E J dR 1 dR 2 F cv (R 1 , R 2 ) 2 ^ IT/pp 

Mi = M| Sl=+ i/2j 1= _3/2;dotl>-+|s2=-l/2J a =+3/2;rfot2> 
~~ JW |si = -l/2,ji=+3/2;dotl>^|s2=+l/2,j2 = -3/2;(ioi2> 

= -£ | dR 1 dR 2 F cv (R 1 , R 2 f- [X ' - ^ • (H) 

Here £ = e 2 d§/e, F av (R u R 2 ) = ^(^^(i?!^^)^^) and d =< >. 
The matrix element < S'lxjX > can also be written as —h/(im )(P cv /E g ), where P cv and 
Eg are the inter-band optical matrix element and the band gap energy of the bulk crystal, 
respectively. 

The matrix element M describes the transfer process with conservation of spin, whereas 
Mi relates to the spin-flip process. The exciton states with J to t = ±2 have no matrix ele- 
ments in our modeli*- The transfer processes with bright excitons (eq. |2J) have the following 
amplitudes: 



M^ x4otl ^b 4ot2 = M + Re[M x ] 
M<p b y ,doti^ii> b y ,dot2 = M Q — .Re [Mi] 
M^b 4otl ^b 4ot2 = -Um[Mi] 
M^ vjdotl ^b )dot2 = iIm[Mi\. (12) 

The matrix elements (fT2|) strongly depend on symmetry. Especially, it is related to the 
off-diagonal transfer processes x <-> y. The off-diagonal amplitudes can be written as 

M i>l,doti^l,dot2 = M lp y! doti^4>b,dot2 = J dR 1 dR 2 F cv (R 1 ,R 2 ) — — 1 ^ -R^\ 5 ~' 

If the double-dot system is symmetric with respect to the inversion operations Rim — > 
— i?i(2), Mi = and the transfer process conserves the linear polarization of excitons, i.e. 
x-exciton in QD1 turns into x-exciton in QD2 and the same rule is applied to y-excitons. 
Therefore spin information can be transferred without losses in the system with spatial- 
inversion symmetry R —R. If the dots are shifted with respect to each other in the xy- 
plane (fig. ^:), spin information in the transfer process becomes partially lost since Mi ^ 0; 
if the shift d is small, Mi cx d x d y . The calculated amplitudes of various transfer processes are 
shown in fig.|U For both QDs, we used the following parameters: L = 2 nm, u^iy) = 
m e = 0.07mo, and rrih = 0.25mo- The above parameters are typical for InAs-based QDs^. 

We note that it is very important to compute the matrix elements beyond the dipole- 
dipole approach since the amplitudes x <-> y vanish within the dipole-dipole approximation. 
Also, the generalized dipole approach (R ^> ai attice ) used in this paper is necessary to obtain 
reliable numbers for all the matrix elements at inter-dot distances R ~ l^ot which are typical 
for experiments. The amplitudes of processes with spin-flip, \if)^,,dotl >— > \ipy,dot2 > and 
\i/)y, dotl >— > \ip^,dot2 >, cannot be obtained within the dipole-dipole approach: Jto[Mi] oc 
R~ 5 for R — > oo. At the same time, the Forster transfer elements with conservation of 
exciton spin has the usual asymptotic behavior at R — > oo: M oc R~ 3 . 

In a normal magnetic field, the Hamiltonian (Q) has an additional term 



HT-7 = M9es z + 9 -fj z )B, (14) 
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FIG. 4: Calculated transfer amplitudes for QD pairs of various parameters and symmetry, a) 
Transfer amplitudes for d = 0. M xx = M yy describe the processes with spin conservation; \M xy \ = 
\Myx\ = 0. (b) Transfer amplitudes at d ^ 0. (c) Calculated transfer amplitudes Mq (spin 
conservation) and Mi (spin- flip) describing transfer processes at high magnetic fields and d = 0. 
In (c), the QDs are strongly asymmetric. 

where B is the normal magnetic field, and g e and g^ z are the g factors. The eigenstates of 
the Hamiltonian in a strong magnetic field are pure states of the total angular momentum: 
\sk = ±1/2, jk = ±3/2 >. Thus, the transfer matrix elements in the limit B — ► oo are 
given by eqs. ^2 The transfer process with conservation of spin is given by the element 
Mo, whereas the spin-flip transfer processes are given be M\. Again, it is important to 
stress the role of symmetry for exciton transfer with spin- flip. If the double-dot system 
is cylindrically symmetric, the spin-flip transfer processes vanish, M\ = 0. In the case of 
asymmetric QDs, Mi ^ 0. Figure |U shows the calculated amplitudes for the case d — 
and strongly asymmetric QDs. One can see that the spin-flip processes become important 
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at small inter-dot distances. Again, the dipole approach would not describe such spin-flip 
effects. 



III. PHONON-ASSISTED COULOMB TRANSFER 

In real QD systems, it is very difficult to find QDs with the same exciton energy. There- 
fore, one should involve acoustic phonons to satisfy the energy conservation requirement 
(fig. EK) • The operator of exciton-phonon interaction reads 



H exc - Ph = j2 xl^^y^^ + *y + ^ qr ' 1 ^ + ct q )} (is) 



— . Here c q is the phonon annihilation operator, r e ^ are the electron (hole) coordinates, 
c ph = 5.6 * 10 5 cm/s is the speed of longitudinal sound, a e ^) = —8.0 eV(1.0 eV) are 
the deformational potentials, and p = 5.3 g/cm? is the crystal mass density. The rate of 
phonon-assisted transfer includes two second-order processes (fig. EJ): 



h _ I < dot2\H exc _ ph \dot2 >< dot2\U C oui\dotl > 



( < dot2\U Co ui\ 


dotl >< dotl\H exc _ p h\ 


dotl > 


-hcphl 







h AE 

5{AE - hc ph \q\), (16) 

where \dotl >= \a;l > |0;2 > and \dot2 >= |0; 1 > \ot';2 > denote the states in which 
an exciton is in QD1 or in QD2, respectively; a = x, y is the spin index of exciton. The 
notation |0; k > means an empty QD. 
Then, the rate (JTfij) is reduced to 

w£ tha _> dot2 , a , = ^ l M ^g^l 2 M F(go)(iV(A ^ + 1)( (17) 

where M c iati,a-^dot2,a' is the inter-dot Coulomb matrix element between the excitonic states 
a and a' given by eqs. IT2l and ITB1 the index a = x,y, q = AE/(hc p h), F(q ) is a function 
given by an integral, and N(AE) is the Bose distribution function at temperature T. 

The rates calculated from eq. fTTI strongly depend on the energy difference AE (fig-EJ). At 
small AE and low temperature T, the rate decreases due to the phonon density of states, 

12 



phonon ^phonon 



Coulomb — I Coulomb 



Dot 1 Dot 2 Dot 1 Dot 2 

FIG. 5: Two contributions to phonon-assisted transfer with different virtual intermediate states. 

2.0 r 




FIG. 6: Calculated rates of phonon-assisted transfer as a function of the resonance AE at T = 4 K . 
The QD pair has inversion symmetry, d = 0. 

whereas at large AE, it becomes small due to the matrix elements of the function e iqr . 
The calculated rate is maximum at AE ~ 2 meV and is about ns' 1 . Since the exciton- 
phonon interaction (fl3j) does not include spin-dependent operators, the spin information 
is not lost in the phonon-emission process. Therefore the spin-selection rules are given by 
the Coulomb matrix elements while the phonon matrix elements conserve the exciton spin 
configuration^. In a symmetric pair of QDs with d = 0, the spin of exciton is conserved. 
However, it is important to note that we have neglected the mixing between heavy and light 
holes; this mixing together with the electron-phonon interaction can result in an additional 
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spin relaxation. In oblate QDs considered here, our approximation is justified since the heavy 
hole-light hole mixing is suppressed due to the large slitting between heavy- and light-hole 
levels in the valence band. 

IV. SPIN-DEPENDENT CROSS-CORRELATION FUNCTIONS 

Similar to the molecular systems^, the coupling between QDs can be seen in the photon 
correlation measurements^. Here we are going to introduce spin-dependent correlation func- 
tions for the case of two coupled QDs. The second-order correlation function is defined as 
g\j =< Ii(t)Ij(t + t) > / < ij(t) >< Ij(t) >, where is the emission intensity of the 

(2) 

i-exciton state. The function g\- is proportional to a number of photon pairs arriving with 
time interval r. 

The nonlinear dynamics of a double-dot system can be quite complex. For simplicity, 
we will consider the limit of weak pumping when the biexciton contribution to the density 
matrix is small. Assuming nonresonant unpolarized excitation of low intensity, we can 
describe the exciton dynamics with a system of linear equations: 

h = T r (n Xtl + n y ,i + n Xt2 + n yj2 ) - 4r P n , 
n X) i = -(r r + T t + T t) s + T s )n X) i + T s n yA + T P n , 
n x ,2 = -(r r + T s )n Xi 2 + T tlS n y ,i + T s n y , 2 + T P n , 

n y>1 = -(r r + r t + r tjS + r s )n x>1 + r s n x>1 + r P n , 

n y ,2 = -(r r + r s )n^ i2 + T^n^i + T s n x ,2 + T P n , (18) 

where no is the "vacuum" exciton state, n a ^ are the numbers of excitons, a = x, y is the 
type of exciton at B — 0, and k — 1, 2 is the QD number. The rate T r describes radiative 
recombination, T t is the energy transfer rate from QD1 to QD2 with conservation of spin, 
T t jS is the inter-dot transfer rate with spin flip, T s is the intra-dot spin-flip rate, and Tp is the 
pumping rate proportional to the light intensity. This simple model resembles the ones used 
in refsj^* 2 ^. For our calculations we choose: T t = 1 ns -1 , T r = 1 Vt, s = T s = 0.1 ns^ 1 , 

and Tp = 0.1 ns^ 1 . The pumping rate Tp = 0.1 ns' 1 corresponds to the regime of low 
intensity (n ~ 1). The small phenomenological spin-flip rates are chosen to take into 
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FIG. 7: Calculated cross-correlation functions for a QD pair with one exciton. The correlation 
function is independent of spin state of excitons. The parameters of relaxation are shown in the 
figure and correspond to a QD pair with R = 5 nm, d = 0, and AE ~ 2 meV. Insert: the energy 
diagram. 

account spin-flip events which are typically slow. The calculated cross-correlation functions 
are shown in figs. I7I81 

(2) (2) 

First we describe one-exciton cross-correlation functions g x {. x 2 and g x ^- y ^ which turn out 
to be spin-independent in our model for the weak pumping regime: g x 2 \. x ^ 2 = 9x\\y,2 = 5l#- 
The time delay r is positive when the photon with energy E exc ^ ti arrives before the photon 
with E excdot2 . At r = 0, emission of the photon E exCjdon projets the system from the state 
\dotl,x > to the "vacuum" state. Therefore, the initial conditions for eqs. EH are set as: 
Hq — 1 and n a ^ k = 0. For r > 0, g{. 2 oc n Xj2 {r) = n^tV), where ftx(V),2 (t) are the solutions of 
eqs. UHlfor the above initial conditions. For r < 0, g{. 2 oc n x> i{r) = %,i(t). The function g\. 2 
is not symmetric with respect to r because of directional exciton transfer from QD1 toward 
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FIG. 8: Upper part: schematics of processes contributing to the correlation function g 



(2) 



. In 



xx,l;x,2 

the first step, a bi-exciton in QD1 emits the x-photon; the second step is inter-dot transfer; the last 

(2) (2) 

step is emission of the x-photon by QD2. Lower part: calculated functions g xa !i x2 an d 9 X xi-y2- 

(2) (2) 

Since inter-dot transfer mostly conserves spin, g x Ji- x2 > 9 X xi-y2- The parameters of relaxation 
are shown in the figure and correspond to a QD pair with R = 5 nm, d = 0, and AE ~ 2 meV. 



QD2. The effective exciton lifetime in QD1 is shorter since an exciton can be transferred 
to QD2. This is reflected as a faster increase of g\.^ at r < 0. In a magnetic field, the 

(2) (2) 

cross-correlation functions can become polarization-dependent (g Xt i- x ^ 9xi- V 2) since the 
resonance conditions can be different for various transfer processes. This can be incorporated 
in the model through appropriate spin-dependent transfer rates F t t j tot -+j> ■ 

Spin transfer processes can be observed using the bi-exciton lines. Even at small pumping, 
weak bi-exciton lines exist . The energy of bi-exciton lines are red-shifted by few me V and can 
be distinguished from the one-exciton lines. In addition, the bi-exciton lines have a quadratic 
power dependence. Consider now the bi-exciton in QD1. It decays in a radiative cascade 
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emitting two photons with the same polarizations (x, x or y, y photon pairs) 24 . If the first 
emitted photon has the x-polarization, the remaining exciton in QD1 has the x-character. 
This exciton can recombine or can be transferred to QD2. Since the exciton spin is mostly 
conserved in the transfer process, we expect that QD2 will strongly radiate x-photons shortly 

(2) (2) 

after emission of the x-photon from QD1. This means that g^xi;x,2 > 9xx,i-y,2> where the 

index xx labels the x-photon emitted by the bi-exciton in QD1. In other words, we use a bi- 

exciton in QD1 as a tool to prepare the x-exciton state at r = 0. Then, the initial conditions 

at r = are: n = 0, n a ^ k = 1 if (a, k) = (x, 1) and otherwise. Figure |H1 demonstrates 

the striking difference between the polarized correlation functions (g xx \i- X) 2 an d i-« 2) i n 

(2) (2) 

the important region r > 0. In this way, by comparing g xxl . x2 and g xx ^ y ^ directional spin 
transfer between QDs can be observed experimentally. 



V. STRONGLY-RESONANT COULOMB TRANSFER 

The convectional Forster mechanism is based on the resonance condition between the 
"donor" and "acceptor"—: the ground-state energy of the donor molecule coincides with the 
energy of an excited state of the acceptor. In the QD system, such a condition can be 
realized if the ground s — s exciton transition in the QD1 has the same energy as the p — p 
transition in QD2 (fig. EJ). Here s and p are the shell indexes in a QD. In this process, the 
s-exciton in QD1 is first transferred to the p-state of QD2; then it relaxes to the ground 
state of the QD2. The transfer rate of this process consists of two contributions: 



Wr o ^ dot 2 = W ph + W dl ^ (19) 

where W ph is the phonon-assisted transfer rate given by eq. EH i n which AE is the energy 
difference between s- and p-excitons in QD1 and QD2, respectively. The rate W dlT describes 
direct resonant transfer between QDs. The latter can be calculated in the spirit of the 
convectional Forster theory as 



2ir 

W Mr = —\M dotl , a ^ do t2,a> \ 2 J{AE), (20) 



h 

1 r en /2 



where J(AE) = j- AE 2/^ r i -n is the normalized effective density of states in the QD2, T en 
is the energy relaxation rate in the QD2, and AE = E excdon — E exc ^ dot2 . To obtain the 
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FIG. 9: Resonant transfer process involving the s-state in QD1 and the p-state in QD2. Relaxation 
in QD2 is phonon-assisted. 

equation ()20j) one should solve the master equation involving the density matrix and assume 
that hY en > \M dotl ^ dot2 \. The latter condition can be easily satisfied because the excited 
p-states are quasi-stationary and the phonon-induced relaxation in QDs is fast usually; the 
typical relaxation times of self-assembled QDs are in the range of 50 pa^. In the opposite 
limit hT en < \M dotl ^ dot2 \, the coupling between the ground s-state of QD1 and the excited 
p-state of QD2 is coherent; it has the character of Rabi oscillations. Then, in the case 
of KT en < \M dot i^ do t2\, the characteristic time for transfer from the s-state of QD1 to the 
s-state of QD2 will be about 1/T en . 

We now calculate the resonant transfer rate for typical parameters of self-assembled QDs. 
Again, the rate strongly depends on the resonance condition (see fig. ITUj). The calculated 
rate demonstrates strong enhancement for AE ~ KT en and also has the structure due to the 
phonon-assisted processes discussed above. 

The transfer rate W res is proportional to \M dot i, a ^ do t2,a' | 2 and therefore the spin-selection 
rules are given by the Coulomb matrix elements. However, the complete transfer process 
contains energy relaxation inside QD2 (fig. EJ). This relaxation can lead to spin flip. Then, 
the efficiency of spin transfer will also depend on the ratio between intra-dot relaxation rates 
with and without spin flip. If energy relaxation inside QD2 involves mostly the heavy-hole 
states, spin-flip relaxation will be weak, since the main contributions to electron-phonon 
scattering are not spin-dependent and, simultaneously, the heavy-hole exciton functions are 
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FIG. 10: Resonant transfer rate corresponding to the process shown in fig. 8. The parameters of 
the QD pair are: R = 5, hw% = 20 meV, hwy = 30 meV, d = 0, and 1/Y en = 40 ps. Temperature 
T = 4 K. The exciton spin is conserved in the transfer process. 

factorized^i. It concerns both the acoustic-phonon interaction (|15J) and the Frohlich scat- 
tering with emission of LO-phonons 20 . To conclude this section, we note that the resonant 
transfer process can involve a localized state in the wetting layer, instead of the p-state in 
QD2. Such a possibility was discussed in refA 



VI. DISCUSSION 



In sections 3 and 5, we discussed incoherent transfer assisted by phonons or involving 
a broadened state in QD2. The phonon-assisted transfer regime between ground states in 
QD1 and QD2 assumes that the energy difference AE is larger than the Coulomb matrix 
element | Md ti~>dot2 1 • Now we briefly consider the case of coherent resonant coupling in the 
regime AE ~ | M^oti^don 1 12 i 13 i 14 . This regime requires fine tuning of energies of QDs; this 
can be done, for example, with magnetic and electric fields. The QDs can be designed from 
different materials and therefore may have different g-factors. By changing a magnetic field, 
one can change AE. A similar principle can be used in the case of an applied electric field; 
if the QDs have different dipole moments, AE can be controlled by the electric field. 
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The calculated Coulomb matrix elements (see fig. @| are in the range of 0.05 meV for 
R = 5 nm. The corresponding time is quite short: At = h/\Mdat\->dot2\ ~ 10 ps. This time is 
much shorter than the typical recombination time of ground-state excitons in QDs, which is 
about 1 ns. Therefore the coherent Coulomb-induced coupling can exist in resonant pairs of 
QDs. In the case of AE = and a QD molecule with d = 0, the one-exciton wave functions 
are given by the linear combinations: 



X 



\dotl, x > \dot2, > ±\dotl, > \dot2, x > 

71 ' 



\dotl,y > \dot2, > ±\dotl, > \dot2,y > 
% = j= • (21) 

The energy splitting within the pairs of states is given by 2|M dotl ^ (io42 |. This energy can 
be regarded as a Rabi frequency. If the exciton is created initially in the QD1, the time 
to transfer the exciton to the neighboring dot would be At = h/\Md a^dot2\ ~ 10 ps for 
the QD pair with R = 5 nm. With increasing R, the transfer times will become longer. 
In the regime of coherent coupling, transfer of spin information occurs coherently, without 
any energy dissipation. The coherent spin-Rabi oscillations between QDs can probably be 
observed with modern optical methods. 

To conclude, we have described the spin-dependent Coulomb interaction in a QD pair. 
Such a coupling suggests the possibility to transfer spin information between individual 
nano-crystals without transfer of charge or mass. The spin-dependent transfer originates 
from the exchange and spin-orbit interactions in semiconductors and strongly depends on 
symmetry and shapes of nano-crystals. If symmetry of a QD pair is high enough, spin 
information can be transferred without losses. If symmetry is broken, spin relaxation in 
the transfer process can become significant. To calculate the transfer rates in the realistic 
model, we use a generalized dipole-dipole approximation which is valid if R 3> ai a tuce- The 
usual dipole-dipole approximation (R ^> Idot) gives too large numbers for the realistic inter- 
dot distances. As a method to observe spin transfer, we consider spin-dependent photon 
correlations in a pair of QDs. 
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